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Abstract. The rate coefficient for radiative and dielectronic recombination of berylliumlike magnesium ions was measured with 
high resolution at the Heidelberg heavy-ion storage ring TSR. In the electron-ion collision energy range 0-207 eV resonances 
due to 2s — > 2p (AN = 0) and 2s — > 31 (AN = 1) core excitations were detected. At low energies below 0.15 eV the 
recombination rate coefficient is dominated by strong 1 s 2 (2s 2p 3 P) 11 resonances with the strongest one occuring at an energy 
of only 21 meV. These resonances decisively influence the Mgix recombination rate coefficient in a low temperature plasma. 
The experimentally derived Mgix dielectronic recombination rate coefficient (±15% systematical uncertainty) is compared 
with the recommendation by Mazzotta et al. (1998, A&AS, 133, 403) and the recent calculations by Gu (2003, ApJ, 590, 1 131) 
and by Colgan et al. (2003, A&A, 412, 597). These results deviate from the experimental rate coefficient by 130%, 82% and 
25%, respectively, at the temperature where the fractional abundance of Mg ix is expected to peak in a photoionized plasma. At 
this temperature a theoretical uncertainty in the Is 2 (2s2p 3 P)7I resonance positions of only 100 meV would translate into an 
uncertainty of the plasma rate coefficient of almost a factor 3. This finding emphasizes that an accurate theoretical calculation 
of the Mg ix recombination rate coefficient from first principles is challenging. 
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1. Introduction 

For the accurate calculation of the ionization equilibrium in as- 
trophysical plasmas, rate coefficients for the population and de- 
population of the various ion charge states have to be known 
precisely. To date most of the required rate coefficients are de- 
rived from theoretical calculations, and, hence, experimental 
benchmarks are required for testing and improving the theo- 
retical methods. This is especially true for recombination in 
photoionized plasmas, which occurs at relatively low plasma 
temperatures of only a few electron volts. At such low tem- 
peratures the recombination rate coefficient, depending on the 
ion under consideration, can strongly be influenced by the ex- 
istence of low-energy dielectronic recombination (DR) reso- 
nances that are difficult to theoretically predict with sufficient 
accuracy. 

Experimentally derived plasma rate coeffici ents were previ- 
ously publishe d for the recombi nation of C iv (jS chippers et alJ 
l200ll) . Ovi teohmetalJl2003l) . Ti v ISchippers et alJll998l) . 
Ni xvm jFode et alJl2003l). Ni xxvi JSchippers et alj|2000l) and 
Fe xvm -Fexxn dSavin et alll997lll999tbo02albll2003l) . Here 
the Mg ix recombination rate coefficient derived from experi- 
mental measurements at a heavy ion storage ring is provided. 
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The approximate temperature ranges where berylliumlike 
magnesium forms in photoionized and in collisionally ionized 
plasma s can be obtained from the work of Kall man & Bau tista 
(200]]) who calculated the fractional abundances of ions in 
plasmas for a variety of physical conditions. For photoionized 
plasmas they find that the fractional Mgix abundance peaks 
at an ionization parameter of log £ = 0.9 corresponding to a 
temperature of about 2.8 eV. The 'photoionized zone' may be 
defined as the range of temperatures where the fractional abun- 
dance of a given ion exceeds 10% of its peak value. For Mg ix 
this corresponds to the temperat ure range 2-13 eV. Using the 
same criterion and the results of Kallman & Bautista ( 200 1|) for 
coronal equilibrium the Mgix 'collisionally ionized zone' is 
estimated to extend over the temperature range 60-170 eV. It 
should be kept in mind that these temperature ranges are only 
indicative. In particular, t hey depend on the accurac y of the 
atomic data base used by iKallman & Bautistal J200ll) and, in 
case of the photoionization zone, on the assumed l/E energy 
dependence of the ionizing radiation. Nevertheless, the above 
given temperature ranges will be used in the discussion below. 

From the present measurements it is found that the Mg ix 
recombination rate coefficient in the photoionization zone is 
decisively influenced by the presence of a strong DR resonance 
at an electron-ion collision energy of 21 meV. Any theoreti- 
cal calculation aiming at an accurate Mg ix rate coefficient will 
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have to predict this resonance's position with an error of less 
than a few meV. Presently available atomic-structure computer 
codes are generally not capable of providing results with such 
an accuracy ab initio. In the standard theoretical approaches the 
situation is often improved by using spe ctroscopic ally obser ved 
target e nergies for the atomic structure JColean et al.l2003tlGul 
2003a). As will be shown below, a 100 meV uncertainty on 
the theoretical resonance position would translate into an un- 
certainty of the Mg ix DR rate coefficient of a factor 2.7 at the 
temperature where the Mg ix abundance is expected to peak in 
photoionization equilibrium. 

2. Experiment 

Various aspects of recombination measurements at the 
heavy-ion storage ring TSR of the Max-Planck-Institut 
fur Kernphysik (MPI -K) in Heid e lberg, G ermany, have 
been described bv [Kilgus et alJ i 19921). Lampert et alJ 
lll996l). IPastuszkaet al] ill 9961) . iMuller & Wolfl Jl997h . and 
ISchippers et alJ J20001 l200ll) . For the present experiment 
a Mgix ion beam with an energy of about 4.2 MeV/u was 
provided by the MPI-K linear accelerator facility and injected 
into the storage ring. There, over a distance of ~1.5 m, the 
ion beam was merged with a collinearly moving, magnetically 
guided electron beam. The electron beam served two purposes. 
First, it acted as a coolant for the ion beam, i. e., the ion 
beam's diameter and momentum spread were reduced by 
collisions w ith the much "colder" electrons (electron cooling, 
IPothl fl990). Second, the electron cooling device (electron 
cooler) was subsequently used as an electron target where 
the ions recombined with the electrons and thereby changed 
their charge state. After adjusting the electron energy to a 
well defined value by tuning the electron cooler's cathode 
voltage appropriately, recombined ions were separated from 
the circulating parent-ion beam in the first bending magnet 
downstream from the electron cooler. Because of their high 
velocity they were efficiently counted with a scintillation 
detector. Absolute recombination rate coefficients as a function 
of relative energy between electrons and ions were derived 
from the measured count rate by i) subtracting the separately 
measured "background" count rate due to electron capture 
from residual gas molecules and by ii) normalizing to the 
number of stored ions and to the electron current which 
were both continuously monitored during the measurement. 
The "merged - beam s technique" was recently reviewed by 
IPhaneuf etallll 19991) . 

For the derivation of absolute rate coefficients from merged 
electron-ion beams measurements the contamination of the ion 
beam by metastable ions is an issue of concern. In a storage- 
ring experiment metastable ions can usually be given enough 
time (up to a few seconds) to decay before data taking is started. 
However, berylliumlike ions with zero nuclear spin like 24 Mg ix 
possess an extremely long-lived metastable Is 2 2s 2p 3 Po state 
which cannot decay to the Is 2 2s 2 So ground state by a one- 
photon transition. DR resonances due to the excitation of the 
3 Pq metastable state were observed in single -pass merged- 
beams experiments with berylliumlike ions bv iBadnell et alJ 
il99lh . Such resonances are not observed in the present ex- 




Electron-ion collision energy (eV) 

Fig. 1. Measured Mg ix merged-beams recombination rate- 
coefficient in the electron-ion collision energy region of DR 
\s 2 2s2pnl resonances attached to Is 2 2s 2 — > \s 2 2s2p 
(AN = 0) core excitations. The dashed line is the calculated 
hydrogenic rate coefficient for radiative recombination (RR). 
Positions of some 1 s 2 2p 2 nl TR resonances are marked by ver- 
tical arrows (see text). 

periment. All measured resonances can be attributed to the ex- 
citation of Mg ix ground-state ions. 

3. Results 

The experimental Mg ix merged-beams recombination rate co- 
efficient is shown in Figures ^ - [5] over different ranges of 
electron-ion collision energy. Figure ^ displays all measured 
recombination resonances due to \s 2 2s 2 — > Is 2 2s2p and 
1 s 2 2s 2 — > 1 s 2 2p 2 (AN = 0) core excitations with the exception 
of the resonances at very low energies below 0.15 eV. These 
were measured separately with higher resolution and are dis- 
played in Figure |2 Finally, Figure [3] shows at higher energies 
the DR spectrum for resonances attached to Is 2 2s 2 — > Is 2 2s 3/ 
(AN =1) core excitations. 

3.1. AN = di- and trielectronic recombination 

Most of the measured AN = photorecombination resonances 
in Figure [2 can be ascribed to Is 2 2s 2 — > Is 2 (2s2p 3 P) 
and Is 2 2s 2 — > ls 2 (2s2p l P) excitations. The corresponding 
Rydberg series of DR resonances converge to their respective 
seri es limits. Their spe ctroscopic values are 17.56 and 33.685 
eV (iMartin et a l. 1999). In order to match the experimental se- 
ries limits - obtained by extrapolating the \s 2 2s2pnl reso- 
nance positions to n — > 00 - to these values a scaling factor 
1.00564 was applied to the energy scale. This factor is within 
the experimental uncertainty. 

Next to the \s 2 2s2pnl DR resonances, further resonances 
are attached to 2s 2 — > 2p 2 double core excitations; they ap- 
pear especially below the ls 2 (2s2p 3 P) series limit. The ra- 
diative decay of these triply excited resonance states to below 
the Mg ix ionization limit completes trielectronic recombina- 
tion (TR). The importance of TR for the photorecombination of 
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Fig. 2. High resolution Mgix merged-beams recombination 
rate-coefficient in the energy range of of the \s 2 {2s2p 3 P)7l 
DR resonances close to zero electron-ion collision energy: 
Experiment (symbols) and RR and DR resonance fit (thick 
full line). Individual contributions to the fit are 6 DR reso- 
nances (thin full lines) and RR (thin dashed line). Data points 
below 3 meV were excluded from the fit. The electron beam 
temperatures resulting from the fit are indicated in the inset 
which shows the same data on a logarithmic energy scale. It 
emphasizes the excess recombination rate coefficient Aa at 
very low energies (see text). The maximum enhancement is 
(Aa + ao)Mo ~ 2.3. Further fit results are given in tabled 



berylliu mlike ions was only recently discovered and discussed 
bv lSchnell et"aT1 J2003h . For berylliumlike CI xiv they find that 
over the plasma temperature range 1-100 eV TR contributes up 
to 20-40% to the total recombination rate coefficient. 

The disentanglement of TR and DR contributions to the 
measured recombination spectrum requires detailed atomic 
structure calculations and is beyond the scope of this paper. 
Therefore, only obvious TR contributions are marked in Figure 
H For a more detailed dis cussio n of TR the reader is referred to 
the work of lSchnell et aDj2003l) . For the present purpose of de- 
riving the experimental Mg ix plasma rate coefficient the ques- 
tion of the origin of individual resonances is not relevant. In any 
theoretical calculation, however, TR has to be accounted for in 
order to arrive at an accurate Mg ix plasma rate coefficient. 

3.2. Recombination at low energies 

Figure |2 shows the experimental merged-beams Mgix recom- 
bination rate coefficient from very low energies up to 0.17 eV. 
In this energy range, the recombination rate coefficient is dom- 
inated by 1 s 2 (2s2p 3 P) 11 DR resonances. For an accurate de- 
termination of the resonance parameters we fitted 6 DR res- 
onance line profiles to the measured spectrum as well as a 
smooth contribution due to radiative recombination (RR). In 
the fit the corresponding cross sections <r(E) were convoluted 
with the experimental electron energy distribution 
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It is characterized by the longitudinal and transversal — with 
respect to the ele ctron beam directio n — temperatures T\\ and 
T ± , respectively dKil gus et alJI 1992b . and £ = (1 - T n /T ± ) 1/2 . 
Furthermore, k B denotes Boltzmann's constant and erf(^) = 
2n~ 1 l 2 C exp -t 2 dt is the error function. 

The convolution of the recombination cross section with 
the above experimental electron energy distribution yields the 
experimental merged-beams rate coefficient 



a MB (E) = ^2E/m e o-(E) f(E, E, T h TJ dE 



(2) 



where E and m e denote the experimental electron-ion collision 
energy and the electron rest mass, respectively. For each DR 
resonance the cross section was taken to be 



r< DR >(£)= Z.h 



F/2 



n E (£ - £ res ) 2 + (F/2) 2 



(3) 



The DR line shape is assumed to be a Lorentzian multiplied 
by a factor E les /E accounting for the 1/E dependence of the 
DR cross section at low energies dSchippers et alJlT998h . For 
E Isfi » T this factor can be neglected. In the limit T — > the 
DR-resonance profile becomes a delta-function and the corre- 
sponding merged-beams rate coefficient can be expressed ana- 
lytically as 



<B fi) (D 



: y/2E K Jm e f(E res , E, T\\, T ± ) 



(4) 



with f(E KS ,E, T\\, T ± ) from Equation^ For T > the merged- 
beams rate coefficient orLg (E) was evaluated numerically. 

For the radiative recombination (RR ) cross section the 
semi-classical hydrogenic formula o f lBethe & Salpeterldl957l) 
was used in slightly modified form, i. e., 



JRR) 



(E) 



2.105 x 10~ 22 cm 2 



^ ] font n 



vis? 



nE(2* s K- 



n 2 E) 
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and the convolution (Eq. |5Jl was performed numerically. In 
equation[5]the letter % denotes the Rydberg constant. The effec- 
tive charge was taken to be Z e ff = 8 and the summation over the 
principal quantum numbers n was carried out from n m ; n = 2 to 
«max = 44. The latter value is determined by field io nization 
in the charge-state-analyzing dipole magnet ( Schi ppers et alJ 
1200 lh . The constants t„ account for partly filled shells. Here 
h - 6/8 and t n = 1 for n > 2 were used. The factors k n cor- 
rect for the deviation of the semi-classical cross sections from 
the more exact quantum mechanical resu lts and were calcu - 
lated following the prescription of And ersen & Bolkol dl990l) . 
They are monotonically increasing for increasing n, starting 
from Ic2 = 0.877 and approaching unity for higher n. 

The resonance parameters that were obtained from the 
fit are listed in Table ^ The reduced y 2 of the fit is 1 .24. 
The electron beam temperatures resulting from the fit are 
k B T\\ = 0.153 ± 0.004 meV and k B T ± = 5.1 ± 0.1 meV. 
These values roughly correspon d to what is e xpected from 
the electron cooler settings (Pastuszka et al. 1996). Data points 
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£ res / [meV] T / [meV] a / [lO" 18 eV cm 2 ] 



20.92(5) 
60.2(1) 
82.5(1) 
99.9(2) 
117.6(9) 
131.0(4) 



4.6(1) 



17.93(7) 
2.23(2) 
1.41(2) 
1.29(2) 
0.30(2) 
0.41(2) 



Table 1. Results of fitting individual DR resonances to the 
measured recombination rate coefficient at low energies (Fig. 
12). Each resonance is characterized by its resonance position 
E les , width r and strength cF. Numbers in brackets denote the 
statistical uncertainties obtained from the fit. If no value for the 
width is given a delta-function-like resonance was assumed. 
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Fig. 3. Measured Mgix merged-beams recombination rate- 
coefficient in the electron-ion collision energy region of 
DR ls 2 2s3lnl' resonances attached to Is 2 2s 2 — > Is 2 2s 31 
(AN =1) core excitations. 

below 3 meV were excluded from the fit. At lower ener- 
gies the measured rate coefficient exceeds the fitted one by 
a factor of up to 2.3 (inset of Figure This recombi- 
nation rate enhancement at very low energies is an inher- 
ent feature of merged-beams experiments at electron coolers 
jGwinner et alJl2000t iHeerlein et al.ll2002l iHorndl et alJl2003l 
IWolf &Gwinnerl20 03). In the present context, this recombina- 
tion rate enhancement can be clearly distinguished from the 
normal photorecombination rate coefficient and is excluded 
from the experimentally derived plasma rate coefficient of 
SectionEH 

3.3. AN = 1 dielectronic recombination 

Resonances attached to Is 2 2s 2 — > Is 2 2s 31 (AN = 1) core ex- 
citations occur in the energy range 70-205 eV (Figure[3j. The 
lowest resonances of this group are the Is 2 2s 31 31' DR reso- 
nances extending to energies of up to about 100 eV. The higher- 
n manifolds of Is 2 2s 31 nl' DR resonances mutually overlap. In 
contrast to AN = DR, the resonance strengths of the AN = 1 
DR resonances rapidly decrease with increasing n. Therefore, 



Fig. 4. Experimental merged-beams recombination rate coeffi- 
cient in the region of high-n Is 2 (2s 2p l P)nl DR resonances 
together with the result of an AUTOSTRUCTURE calcula- 
tion scaled by a factor 0.7. The dashed line is the theoreti- 
cal result with account for the experimental field ionization 
of high-» Rydberg sta tes using the field-ionization model of 
Schip pers et aH feOOll) . The full curve is the theoretical result 
including the full DR resonance strength up to n = 1000. 
The inset shows the same curves on a different scale for an 
overview. 

the resonance strength does not pile up at the series limits. 
In Figure |3 two series limits can be discerned at 202.25 and 
205. 14 eV. In order to match these valu es from the NIST atomic 
spectra data base (Mar tin etall fl999) for the \s 2 2s3d l Dnl 
and the Is 2 2s 3d 3 Dnl series, respectively, the experimental 
energy scale was multiplied by a factor 1.0091. This rescaling 
is within the experimental uncertainty. 

3.4. Plasma recombination rate coefficient 

Formally, the recombination rate-coefficient in a plasma is ob- 
tained by using an isotropic Maxwellian distribution function 
in equation|2l i. e., 

a p iasmafe7')= .- \ ■ f (r(E)E exp ) dE. (6) 

For plasma temperatures k B T » k B T ± a factor cr(E) y/2E/m e 
in the integral may be replaced by the experimental merged- 
beams recombination rate coefficient. This approach can safely 
be used for all recombination resonances with E [es » k B Tj_ . 
Here it was applied to all recombination resonances above 0.2 
eV that are shown in Figures ^ an d El I n order to exclude any 
effect of the finite experimental energy spread on the experi- 
mentally derived plasma rate coefficient, Equation |S] was used 
with the resonance cross section from Equation|5]for the reso- 
nances at lower energies. The contribution of the recombination 
resonances below 0.2 eV (Figure |5Jl to the plasma recombina- 
tion rate-coefficient was evaluated by using the fitted resonance 
parameters from Tabled 

For the derivation of the experimental Mg ix recombination 
rate coefficient one has to be aware of the fact that loosely 
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i a I [cm 3 s" 1 K 3/2 ] 



E, I [eV] 



1 


3.5960 x 10 


2 


2.0071 x 10 6 


2 


2.0109 x 10" 


-2 


3.7011 x 10 s 


3 


2.6244 x 10" 


-3 


1.2038 x 10 5 


4 


3.6203 x 10" 


-4 


4.5163 x 10 4 


5 


4.1121 x 10" 


-5 


9.2881 x 10 2 


6 


2.6915 x 10" 


-5 


2.2683 x 10 2 


7 


6.7347 x 10" 


-7 


3.7786 x 10 1 



Table 2. Parameters for the fit of equation to the experi- 
mentally derived Mg ix DR+TR rate coefficient in a plasma. In 
the temperature ranges 25-350 K and 350-9 x 10 8 K the fit 
deviates less than 2% and 1%, respectively, from the experi- 
mentally derived result. The systematic uncertainty of the ex- 
perimentall y derived absolute recombination rate coefficient is 
+ 15% (Lampe rtetalJl996l) . 



bound high-n Rydberg states are easily field ionized in our 
experimental setu p. This effect was modelled in considerable 
detail by Schipper s et alJ J200ll) . Generally, it is more signif- 
icant for lower charged ions. As shown in Figure |4] a theo- 
retical calculation of the Mg ix PR rate co efficient using the 
AUTOSTRUCTURE code of lBadnelll Jl986l) can only be rec- 
onciled with the experimental data if field-ionization is taken 
into account. 

In the present case the effect of field-ionization on the mea- 
sured DR spectrum can be approximately described as a cutoff 
of high-n Rydberg resonances with n > 44. In order to account 
for the unmeasured DR resonance strength the DR spectrum for 
n < 44 is extrapolated to higher n by scaling the result of the 
AUTOSTRUCTURE calculation to the experimental spectrum 
in the energy range 28.0-32.9 eV (scaling factor 0.7, Figure^} 
and by replacing the measured merged-beams rate coefficient 
by the scaled calculated value in the energy range 32.9-34. 1 eV. 

In principle, such a correction also has to be made for all 
other series of Rydberg resonances. For AN = DR such other 
series are relatively weak (Figure [0 and their cutoff by field- 
ionization does not lead to a significant error when consider- 
ing the overall +15% experimental uncerta inty of the absolute 
rate coefficient scale dLampert et alj fl996). Regarding AN - 1 
DR, the resonance strength, as already mentioned, decreases 
strongly with increasing n, leading again to a negligible effect 
of the experimental field-ionization on the plasma rate coeffi- 
cient. 

Our experimentally derived Mg ix DR+TR recombination 
rate coefficient in a plasma — including the extrapolation of the 
1 s 2 (2s 2p 1 P) nl resonance strength to n = 1000 (cf. Figure|4} 

— is shown in Figure [5] as a thick full line. For the convenient 
use of our result in plasma modelling codes the following func- 
tional dependence — customarily used for DR rate coefficients 

— was fitted to our experimentally derived curve: 
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Fig. 5. Comparison of the present experimentally derived 
Mgix D R rate coefficien t (thick full line) with the recommen- 
dation of lMazzotta et all dl998[ dash-dot-dotted line) and the 
theore tical results o flGuH2003al dashed line) and lColgan et alJ 
(2003, dash-dotted line). The error bars denote the ±15% ex- 
perimental uncertainty of the absolute rate coefficient. The thin 
full line is the RR rate coefficient calculated (eq. |5Jl with the 
RR cross section from equation|5]with n max = 1000 and the re- 
maining RR parameters from section [3~2l The dotted line is the 
theoretical RR rate coefficient of lGul(l2003bl) . The temperature 
ranges (see section [0 where Mgix exists in photoionized and 
in collisionally ionized plasmas, respectively, are highlighted. 



In the fit to the low-energy experimental merged-beams re- 
combination rate coefficient (Figure |2J RR was also included 
(with T\\ and T ± as the only free fit parameters and all other 
parameters kept fixed). The same RR cross section was used to 
draw the RR contribution in Figure^ This shows that Equation 
|5]for the RR cross section is consistent with our experimental 
data. Therefore, Equation |5] in conjunction with Equation [6] 
can also be used to derive the experimental Mgix RR rate co- 
efficient in a plasma. 

To this end, the RR cross section parameters were taken 
from section l3~2l In order to make up for field-ionization effects 
the plasma RR rate coefficient in Figure [5] (thin full line) was 
calculated with n max = 1000 instead of n max = 44 as in Figure 
12 For convenient further use our Mg ix RR r ate coefficient in a 
plasma was fitted by a formula introduced bv lVerner & Ferlandl 
(1993): 



o< (RR) - A 



1+J^ 
To 



l-b 



Ti 



l+b 



(8) 



During the fitting procedure the coefficients c, and Ej were var- 
ied. The fit results are listed in Table[2] 



The parameter values that were obtained from the fit are A = 
5.7414 x lO" 10 cm 3 s"\ b = 0.67692, T = 206.15 K, and 
Ti = 7. 6200xl0 6 K. Over the temperature ranges l-2.9xl0 6 K 
and 2.9 x 10 6 - 10 8 K the fit is accurate to more than 2% and 
3%, respectively. 

In Figure|5]the present experimentally derived RR rate co- 
effici ent is compared with the the recent theoretical result of 
iGvJ d2003bl dotted line). The difference is less than 17% over 
the temperature range where Mgix exists in a photoionized 
plasma. 
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The total (unified) Mg ix recombination rate coefficient is 
readily obtained as the sum of our DR+TR and RR rate co- 
efficients. It should be noted that this result does not depend 
on the existence of possible quantum mechanical interferences 
between RR and DR, since interference has already been ne- 
glected during the decomposition of the measured spectrum 
into DR and RR (Figure |2j. In principle, interference between 
RR and DR can lead to asymmetric recom bination resonance 
line shapes (see e. g. ISchippers et alJE()02h . No evidence for 
such effects is found in our experimental Mgix recombina- 
tion spectra. The systematic uncertainty of the experimentally 
derived total Mgix recombination rate coefficient is ±15% 
jLampertet alll996l) . 

4. Discussion and Conclusions 

In Figure [5] the present experimentally derived DR+TR 
rate coefficien t is co mpared with the recommendation of 
iMazzotta et alJ dl99Rh and with the recent calculations of iGvJ 
J2003al) and lColgan et aT]d2003h. T h e reco mmended Mg ix DR 
rate coefficient of iMazzottaetal] lll998h deviates by up to 
130% from the experimental result in the temperature range 
2-13 eV where Mgix exists in photoionized plasmas. This 
rather large discrepancy questions the usefulness of the rec- 
ommended Mg ix DR rate coefficient for the modeling of pho- 
toionized plasmas. In the temperature range 60-170 eV (colli- 
sionally i£nized_zon£Uherecommended Mg ix DR rate coeffi- 
cient of Mazz otta et alJ dl998l) agrees within 11%, i. e., within 
the 15% experimental uncertainty, with the experimentally de- 
rived one. 

In this t emperature range the theoretical rate coefficient of 
IGuI d2003ah shows very good agreement with the experimen- 
tally derived rate coefficient. The difference between the two 
curves in this range is less than 7%. At lower temperatures, 
however, the discrepancy is up to 99% at 2 eV. At these low 
temperatures where Mgix exists in photoionized plasmas, the 
theoretical result o f Colgan et al. (2003) agrees better with the 
experimental rate coefficient. The deviation is less than 35% 
for ksT > 2 eV, less than 30% for kgT > 15 eV, and less than 
15% for k B T > 100 eV. At the temperature 2.8 eV where the 
fractional abundan ce of Mg ix is expected t o peak in a pho- 
toionize d plasma ( Kall man & Bautistall200ll) the results of lGvJ 
( 2003_a|) and lColgan et alJ ( 120031) deviate from the experimental 
rate coefficient by 82% and 25%, respectively. 

The rather large l ow-t emperature deviation of the theoreti- 
cal rate coefficient of Gu (2003a) from the experimentally de- 
rived curve is most probably due to an inaccurate calculation of 
the DR resonance energies below 0.15 eV (Figure|2j. At these 
very low energies the DR rate coefficient is very sensitive to 
slight variations of resonance positions. This is highlighted in 
Figure [6] that displays the effect of hypothetical shifts of the 
low-energy resonance positions by ±50 meV and ±100 meV 
with respect to the tabulated values (Tabled on the Mg ix DR 
rate coefficient. Obviously, a ±100 meV uncertainty of the low- 
energy Is 2 (2s 2/? 3 P) 11 resonance positions translates into an 
uncertainty of a factor 2.7 of the plasma rate coefficient in 
the temperature range where the fractional abundance of Mg ix 
peaks in a photoionized gas. 
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Fig. 6. Impact of a shift of the low energy Is 2 (2s 2p 3 P)ll 
resonances displayed in Figure |2] on the Mg ix DR plasma 
rate-coefficient. The full curve is the experimental result, the 
long-dashed curves correspond to resonance shifts of AE KS = 
±50 meV, and the short-dashed curves to a shift of AE res = 
±100 meV. Thereby, the lower curves correspond to resonance 
shifts to lower energies. The vertical dash-dotted line marks the 
temperature 2.8 eV where th e fractional abundance of Mgix 
peaks in a photoionized gas jKallman & Bautistall200ll) . The 
horizontal dash-dotted lines indicate the factor 2.7 uncertainty 
of the plasma rate coefficient at this temperature that would be 
introduced by a ±100 meV uncertainty of the low energy reso- 
nance positions. 

These findings demonstrate that an accurate calculation of 
DR rates for ions that exhibit DR resonances close to zero en- 
ergy is challenging. In case of more complex ions this task is 
certainly beyond the capabilities of present-day atomic struc- 
ture codes as is exemplified by a recent combined theoretical 
and e xperimental DR study of argon-like Sc 3+ ( Schi ppers et alJ 
2002). Recombination experiments at heavy-ion storage rings 
are certainly required for guiding the future development of the 
theoretical methods, especially in the case of low temperature 
DR rate coefficients for complex ions. 
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